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Dangling bond wires on Si(OOl) are prototypical one dimensional wires, which are expected to 
show polaronic and solitonic effects. We present electronic structure calculations, using the tight 
binding model, of solitons in dangling-bond wires, and demonstrate that these defects are stable in 
even-length wires, although approximately 0.1 eV higher in energy than a perfect wire. We also note 
that in contrast to conjugated polymer systems, there are two types of soliton and that the type of 
soliton has strong effects on the energetics of the bandgap edges, with formation of intra-gap states 
between 0.1 eV and 0.2 eV from the band edges. These intra-gap states are localised on the atoms 
comprising the soliton. 

PACS numbers: 68.65. -hk; 73.20.-r; 71.15.Nc; 71.38. -hk 



I. INTRODUCTION 

There are many motivations to understand the 
transport properties of materials in the extreme one- 
dimensional limit. Some are technological: the logi- 
cal conclusion to the historic reduction in size of elec- 
tronic components would be active device elements, and 
passive connections between them, that are of atomic 
scale. There has been much recent interest in structures 
that might act as atomic-scale wiresi or as atomic- or 
molecular-scale switches. Other reasons relate to fun- 
damental physics: transport in one dimension is very 
different to higher dimensions, because the coupling of 
electrons to other excitations, both of the lattice and of 
the electronic system, is strong. This can lead to insta- 
bilities such as the formation of a Luttinger liquid^ or one 
of a Peierl92i (infinite wire) or Jahn-Telle»^ (finite wire) 
distortion, with correspondingly strong modifications to 
transport properties. 

One striking example of this type of behaviour is in 
conductive polymer systems^. These exhibit an alterna- 
tion of double and single bonds which is at least partly 
of the classic Peierls type (driven by electron-lattice cou- 
pling), although electron correlation effects may also be 
important^. Carriers introduced into these systems are 
localised, probably partly by disorder but also to a signif- 
icant degree by self-trapping by the lattice^. There are 
two main classes of charged defect; one is the polaron, 
which is similar in principle to excitations in higher- 
dimensional systems. Here the carrier is surrounded by 
a region in which the atomic distortions serve to lower 
its own energy (in this case, involving a reduction in the 
double-single bond alternation). The second type is the 
soliton; this can occur only in 'degenerate' systems, such 
as irans-polyacetylene (t-PA) , where there are two equiv- 
alent ground states of an infinitely long chain, in which 
the double and single bonds are interchanged. The soli- 
ton involves a mid-gap state associated with a 'domain 
wall' between regions of opposite bond alternation. Al- 
though localized, both polarons and solitons are highly 
mobile, and much is known about their effect on trans- 
port in both the coherent and incoherent limit a^i^i^ . 



In this paper, we examine an alternative pseudo-one- 
dimensional system based on the "dangling-bond wire". 
This system is formed on a hydrogen terminated silicon 
(001) surface via STM-induced selective desorption of hy- 
drogen along the edge of a dimer row^ft. The atoms in 
such wires have been shown (both theoretically and ex- 
perimentally) to undergo a Peierls/ Jahn- Teller type dis- 
tortion, giving rise to an alternating pattern of atomic 
positions with respect to those of the passivated sur- 
face and producing the electronic structure of a one- 
dimensional narrow-gap semiconductopifi. In contrast to 
the conducting polymers, this distortion occurs predom- 
inantly perpendicular to the axis of the 'wire'; however, 
as in the conductive polymers, the resulting structure is 
not expected to be "static". We recently predictedii, on 
the basis of tight-binding calculations, that electrons or 
holes introduced into such wires will produce self-trapped 
"small polaron" defects. Despite their localisation, we 
have also shown that these defects are remarkably mobile 
near room temperature and aboveiA, raising the interest- 
ing possibility of the transport of charge through such 
devices. 

It is natural to ask whether there might be solitonic, 
as well as polaronic, defects in dangling-bond wires (as 
in t-PA); after all, these structures are degenerate (the 
total energy is invariant under exchanging the "up" and 
"down" atoms) . This question is intimately related to the 
effect of different types of boundary condition; in the case 
of t-PA, double bonds are "anchored" at the ends of the 
chain, so that the ground state of a neutral chain with 
odd length contains a bond-alternation defect (soliton) 
at the centre. The behaviour in t-PA is simple because 
the energy scales for the binding of double bonds to the 
ends of the chain are much larger than any involved in 
the formation of defects along it. 

In this paper we explicitly focus on the behaviour of 
dangling-bond wires of finite length, with a view to study- 
ing the formation of solitons and the end effects. We find 
that the situation is more complicated than in t-PA, be- 
cause the energy scales of the end effects and the defects 
are comparable to one another. This means that there 
is much more complex behaviour when deliberate defects 
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are introduced into the alternating pattern. Accordingly 
we pay particular attention to the stability of various 
possible defects and their contributions to the electronic 
structure of the wire. In light of the large number of 
possible initial configurations and necessarily large sys- 
tem size (>400 atoms), we chose to use a semi-empirical 
tight-binding approach'? that was fast enough to allow 
a comprehensive study of possible defects and had ac- 
curate parametrisations available for the system under 
consideration. 



II. TECHNICAL DETAILS 
A. Basic configuration & relaxation parameters 

A silicon slab was constructed out of 6 layers, paral- 
lel to the (001) plane. Each layer had a 4 x 12 array 
of silicon atoms positioned at the appropriate positions 
for a perfect "diamond" lattice, with the "bottom" layer 
terminated with hydrogen atoms. This layer and the as- 
sociated hydrogen atoms were not allowed to move during 
relaxation and represented the bulk. The unterminated 
upper surface was allowed to reconstruct to form a pair 
of dimer rows. This (001) surface was then terminated 
with hydrogen atoms to produce a fully terminated slab 
that served as a basis for all the simulations. This unit 
cell contained 432 atoms. 

Each slab was separated from its neighbours by a vac- 
uum gap of 15 A once periodic boundary conditions were 
applied. For the maximum length considered of eight 
dangling bonds, each wire was isolated by four hydrogen 
terminated atoms from its virtual neighbours along the 
wire direction, while shorter wires had proportionately 
greater numbers of hydrogen terminated atoms at each 
end. All wires were isolated by a completely terminated 
dimer row from adjacent wires once periodic boundaries 
were considered as this has been found to provide ade- 
quate isolationiA. 

Wires were created by removing the required number 
of terminating hydrogen atoms from the upper surface, 
along one side of a dimer row. The reference point for 
distortions was the average position of hydrogen termi- 
nated silicon atoms in the passivated surface, labelled as 
the baseline or "level" position. Displacements normal 
to the surface in the direction of the bulk were "down" 
and away from the bulk were "up". Before relaxation, up 
atoms were displaced by an additional 0.3 A compared to 
terminated atoms of the normal spacing and down atoms 
by 0.4 A less than terminated atoms. 

Previous worlsi^ had examined the effect of slightly 
perturbing the starting positions to determine whether 
the final results were stable with respect to small dis- 
tortions and found this was the case. Accordingly only 
one magnitude of starting displacement was considered 
in this study. 

All relaxations and electronic structure calculations 
were performed using the Oxford Order-N (OXON) tight- 



Mesh Size 


3 atom chain 


8 atom chain 


1x1x1 


-1963.7827 


-1948.6625 


1x2x1 


-1963.7828 


-1948.6625 


2x1x1 


-1964.1376 


-1949.0197 


2x2x1 


-1964.1377 


-1949.0197 


2x4x1 


— 1964.1377 




4x1x1 


-1964.1364 


-1949.0183 


4x2x1 


-1964.1365 




4x4x1 


-1964.1365 




4x4x4 


-1964.1365 





Table I: Total energies (eV) calculated for relaxed perfectly 
ordered 3 and 8 atom long chains at K as a function of 
fc-point mesh size. 



binding packagei^. Structural relaxations were per- 
formed at a system temperature of K until the maxi- 
mum force per atom reached the limit of 0.01 eV/A. The 
majority of simulations required in the order of 150 iter- 
ations to relax. The Hamiltonian was solved using exact 
diagonalisation and periodic boundary conditions were 
applied to the simulation. Spin-polarisation was not in- 
cluded in the calculations. The silicon parametrisation 
used was that of Bowler et al^, which was specifically 
fitted for the Si(OOl) surface and hydrogen on that sur- 
face. It is worth noting that this parameterisation did 
not fit to the conduction band, so that absolute values 
should not be trusted, though generic behaviour (such as 
shifts in levels) are likely to be correct. 

Convergence of the total energy was tested with re- 
spect to the fc-point mesh. An initial search was per- 
formed with a 3 atom wire to minimise calculation time, 
and selected results replicated with an 8 atom wire. Ta- 
ble ffl shows the results of the convergence calculations. 
Results were adequately converged (<0.01 eV) using a 
Monkhorst-Packii 2x1x1 grid, and all simulations pre- 
sented below used such a mesh. 



B. Configurations considered 

Dangling bond wires 4, 6 and 8 atoms long were con- 
sidered. Some simulations of odd length wires were per- 
formed, but these were not pursued for reasons discussed 
below in Section ITlI Al Each wire was considered initially 
with perfect structure, i.e. the initial starting positions 
were alternately up then down (or vice versa). Defects 
were then introduced consisting of two adjacent up atoms 
or two adjacent down atoms at all possible positions in 
the wire with the remaining atoms alternating appropri- 
ately. 



Figure 1: A DB wire of length eight, with a DD soliton in the middle (marked with darker atoms). The image shows only the 
top four layers of silicon, and limited terminating hydrogens (one on the left and two on the right). 



III. RESULTS 
A. End Effects 

We have previously tested the use of tight binding for 
infinite dangling bond wiresU^, obtaining good agree- 
ment with experiment and ab-initio results. However, 
when we came to examine finite wires with an odd num- 
ber of sites, we found an anomaly: one site on the wire 
became "level" (i.e. at the same height as if it were hydro- 
genated) and the rest of the wire behaved as a wire with 
an even number of sites, in contrast to ab-initio results. 
The local charge on the level atom was zero, again in con- 
trast to the rest of the wire (where down atoms showed a 
charge deficit and the up atoms showed a charge excess) . 
In this section, we explore the reasons for this behaviour 
and investigate whether tight binding is a viable tech- 
nique for such systems. 

The shortest, and simplest, finite wires of odd and 
even length for which there is experimental data are 
those with three and four sites respectiveljiifi (the length 
two wire might be too short, as it effectively consists of 
only end atoms). We have already performed extensive 
investigations of this system using spin-polarised den- 
sity functional theory, and found good agreement with 
experimenti^. In particular, we found that the length 
three wire has a rather small displacement pattern (with 
the two end atoms 0.11 A higher than the central atom) 
in close agreement with experiment (where displacements 
of 0.15 A are measured), and that the length four wire 
can form alternating up/down or down/up patterns, or 
an up/down/down/up pattern. 

Our tight binding simulations of even length wires were 
in good agreement both with the experimental results 
and the DFT simulations. In order to better understand 
the peculiar behaviour for the odd length wires, we used 
local charge neutrality (LCN)i2, (where the on-site en- 
ergies of individual atoms are adjusted until their net 
charge is zero) to remove any charge transfer effects in the 
tight-binding simulations (it is worth noting that it is still 
perfectly possible to obtain buckled dimers on the clean 
Si(OOl) surface with a LCN condition — we will consider 
the effects on an even length wire below). The length 
three wire became almost flat under this condition — 
though the two end atoms were 0.06 A higher than the 



central atom, in qualitative agreement with experiment 
and DFT modelling. We also modelled the length four 
wire with LCN, and found little change — the alternating 
up/down pattern was still observed. However, the LCN 
simulations were extremely hard to converge, and did not 
show particularly improved results. We conclude that, 
for odd length wires, there is spurious charge transfer 
that makes modelling of these systems potentially inac- 
curate in tight binding, and we will not consider them 
furthei"^. However, for even length wires there is no 
problem and we will continue to model these using tight 
binding with confidence. 



B. Solitons in even length wires 

Simulations were performed for chain lengths of 4, 6 
and 8 dangling-bond systems. Initial configurations were 
chosen so as to include all possible permutations of sin- 
gle position ordering defects, for both up and down types. 
All relaxations successfully converged to the force limit. 
The results are shown in detail in Table Hll but in sum- 
mary both UU and DD soliton defects could be formed 
in the relaxed wire. Compared to a perfectly ordered 
wire, those containing solitons are approximately 0.1 eV 
less stable. No soliton defects formed on the end atoms 
and there was strong periodicity in possible locations of 
the defects, as seen in the table. 

As well as showing the final, relaxed configuration of 
the wires (with up atoms notated "u", down atoms "d" and 
atoms associated with a soliton in bold capitals) and var- 
ious energies, Ta.ble UTI shows the initial configuration for 
each system. It is important to note that several different 
initial configurations led to the same final configuration 
(for reasons discussed below). 

The effects of the ends of the wire are rather small 
(only a few meV) but seem to exert a small attractive 
force on the solitons. More interesting is the apparent 
instability of solitons of both kinds exactly at the ends 
of the wire (and of certain configurations with the soli- 
tons in the middle of the wire). For instance, the third 
and fourth lines of table lul show that the initial config- 
urations DDud and UUdu relaxed back to perfect wires, 
while the eighth line shows that duDD relaxed back to 
dUUd. This pattern can be easily understood in terms 
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Table III: Values of the bandgap in even-length wires con- 
taining various defects. VB stands for Valence Band, CB for 
Conduction Band, all energies in eV. 
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Figure 2: Energies (eV) of states adjacent to the band gap 
for a) perfectly ordered 8 atom wire, b) 8 atom wire with UU 
defect, c) 8 atom wire containing DD defect. 



Table II: Final structure, energy change, valence and conduc- 
tion band offsets and starting structure for soliton-containing 
even-length wires. Edifi represents difference in total en- 
ergy relative to a perfect wire. All energy differences are in 
electron-volts. Soliton defects are shown in bold capitals. VB 
stands for valence band, CB for conduction band. 



of charge balance: an up atom is associated with an ex- 
cess of charge, while a down atom is associated with a 
deficit of charge. The configurations which were unstable 
(which includes all systems with a single soliton at the 
end of the wire) had different numbers of up and down 
atoms, which would lead to charge imbalance, so they 
changed to a stable configuration with equal numbers of 
up and down atoms, either removing the soliton entirely 
or changing it. 

Looking at the electronic structure of relaxed wires 
containing solitons, the presence of a soliton reduced the 
bandgap. The general pattern is that a DD type defect 
lowered the bottom of the conduction band, whilst an 
UU type defect raised the top of the valence band. The 
values for different length wires are shown in table IIIll If 



the soliton was of the DD type, the bottom of the con- 
duction band was lowered by 0.1 eV, whilst the top of the 
valence band was raised by 0.22 eV by the presence of an 
UU defect. Both these effects were due to the introduc- 
tion of a intra-gap state, as illustrated in Figure El We 
recall that the system was not parameterised for the con- 
duction band and accordingly absolute energies relating 
to it should be treated with caution. 

A soliton is generally associated with a mid-gap state, 
which is seen here in the change of the valence and con- 
duction bands. These effects can be understood in terms 
of the hybridisation of the silicon making up the surface: 
in ideal, bulk positions the atoms are sp'^ hybridised, but 
reconstruction into dimers at the surface pulls the atoms 
away from this state. When Jahn- Teller distorted, either 
during buckling of dimers or forming a finite length DB 
wire as here, the atoms displaced up move closer to sp'^ 
hybridisation and gain charge (tending to a lone pair) 
while the atoms displaced down move closer to sp^ hy- 
bridisation with a lone p-orbital and lose charge (tending 
to an empty dangling bond). Thus the up atoms are as- 
sociated with filled states at the top of the valence band. 
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a) VB dudUUdud 

b) CB uduDDudu - 



Atom position (Angstroms) 

Figure 3: Effects of soliton formation on the electronic struc- 
ture of finite DB wires. We plot contributions to the charge 
density for specific bands from atoms most strongly affected 
as a function of distance along the wire; a) highest occupied 
electronic state for a central UU soliton, b) lowest unoccupied 
electronic state for a central DD soliton. 



while the down atoms are associated with empty states 
at the bottom of the conduction band. The up and down 
displacements affect the underlying substrate (in partic- 
ular the second and third layer atoms) in such a way that 
alternating up and down displacements (either along or 
across a dimer row) are energetically favourable. When 
this is interfered with, as in a soliton, the extent of the 
relaxation towards sp^ or sp^+p is reduced, and the band 
edges are affected. This is in marked contrast to conju- 
gated polymers, as will be discussed below. 

The solitons are broadly locaHsed as expected, which 
is shown in Figure El This is similar to the polaron seen 
in the same syste m^^'^^ and indicates that the soliton is 
similarly weakly coupled to the bulk. The overall be- 
haviour and form of the soliton in the DB wire is very 
different to the conjugated polymer system. In particu- 
lar, the soliton is associated with a pair of atoms and has 
two forms, while in conjugated polymers it is associated 
with the bond lengths near a single atom, and has only 
one form. 



IV. DISCUSSION AND CONCLUSIONS 

Although the dangling-bond system is too reactive to 
be a component of nanoscale circuitry, it does serve both 
as a potential model for the study of one-dimensional 
conductors and as a stepping stone to realisable useful 
systems, for example as a template for the deposition of 
metal atoms in linear features. As such it is important 
that systems can be modelled accurately and quickly. 
In this paper, we have demonstrated that tight-binding 
simulations are suitable for even-length wires and used 
such techniques to explore the stability of ordering de- 
fects in the Jahn- Teller alternating displacement seen in 
perfect dangling bond wires. We note that modelling 
of odd-length wires is complicated by the development 
of "level" atoms in the relaxed system. This appears 
to be related to difficulties in adequately modelling un- 
paired electrons using the tight-binding formalism; we 
have modelled these systems using spin-polarised DFT, 
reported elsewhere^^. 

In the even length system, we demonstrate that "soli- 
ton" type defects, where alternation between up and 
down displacements is interrupted, are only ^ 0.1 eV 
less stable than the perfect wire, but are associated with 
strong effects on the edges on the valence and conduction 
bands, as shown in Figure|2l as a consequence of changes 
in hybridisation state. The presence of an "UU" defect 
leads to the formation of an isolated state 0.22 eV above 
the original valence band. A "DD" defect is associated 
with a state 0.1 eV below the original conduction band, 
although the conduction band energies should be treated 
with caution due to the nature of the parametrisation. 
We also demonstrate that these states contain large con- 
tributions from the atoms associated with the defect i.e 
localisation. 
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